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Abstract 

In confined systems near a continuous phase transition the long-ranged fluctuations of the corre- 
sponding order parameter are subject to boundary conditions. These constraints result in so-called 
critical Casimir forces acting as effective forces on the confining surfaces. For systems belonging 
to the Ising bulk universality class corresponding to a scalar order parameter the critical Casimir 
force is studied for the film geometry in the crossover regime characterized by different surface 
fields at the two surfaces. The scaling function of the critical Casimir force is calculated within 
mean field theory. Within our approach, the scaling functions of the critical Casimir force and of 
the order parameter profile for finite surface fields can be mapped by rescaling, except for a narrow 
crossover regime, onto the corresponding scaling function of the so-called normal fixed point of 
strong surface fields. In the crossover regime, the critical Casimir force as function of temperature 
exhibits more than one extremum and for certain ranges of surface field strengths it changes sign 
twice upon varying temperature. Monte Carlo simulation data obtained for a three-dimensional 
Ising film show similar trends. The sign of the critical Casimir force can be inferred from the 
comparison of the order parameter profiles in the film and in the semi-infinite geometry. 
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I. INTRODUCTION 



Finite-size contributions to the free energy of a fluid confined between two planar surfaces, 
separated by a distance L, give rise to an effective force per unit area between the surfaces, 
or an excess pressure. This so-called solvation force fsoiv depends on L, the thermodynamic 
state of the bulk fluid, the fluid-fluid interactions, and the two substrates potentials jl|. At 
the bulk critical point of the fluid the solvation force acquires a universal contribution which 
is long-ranged in L. This fluctuation induced effective force is called the critical Casimir 



force 



The critical Casimir effect is a subject of considerable theoretical and experimental in- 
terest, involving experiments for wetting films near critical end points and for colloidal 



systems in the presence of a critical solvent 



-|l2[|. The sensitive temperature dependence 



of the critical Casimir force can be exploited in the latter systems in order to control the 
collective behavior of colloidal particles, such as their aggregation behavior, which opens 
up application perspectives in many areas of material science. According to the accumu- 
lated knowledge, the sign of the critical Casimir force can be selected by suitable surface 
treatments j3,Il3|, 3- Recently, a continuous tuning has been achieved experimentally for a 
colloidal particle in a critical solvent and near a substrate with a gradient in its preferential 



adsorption properties for the two species forming the binary liquid mixture as a solvent |15 |. 
It is very encouraging that the adsorption preference of a substrate can be changed contin- 
uously between strong adsorption of one species to strong adsorption of the other species of 
a binary liquid mixture by tuning the chemical composition of a monomolecular overlayer 
only, without altering the bulk material of the confining substrates. As will be discussed 
later, in the present context this amounts to continuously tune a surface field hi, which 
expresses this preference and breaks the symmetry of the order parameter, between -|-oo 
and — oo. 

Such a tunability of critical Casimir forces towards repulsion might be relevant for micro- 
and nano-electromechanical systems in order to prevent stiction due to the omnipresent 
attractive quantum mechanical Casimir forces 16| - for exceptions see Ref. 17|]. In order 
to achieve repulsive quantum Casimir forces, rather complex systems have been considered 



but they are not yet experimentally established 



181]. Here we show theoretically that, for 



suitably prepared system parameters, switching the sign of the critical Casimir force can be 



2 



achieved not only by varying the surface fields but also via minute temperature changes. 
This occurs if the adsorption properties of the confining surfaces differ significantly, for 
example, if each surface attracts a different component of the binary liquid mixture but one 
does so weakly and the other strongly. 

To be specific, we investigate theoretically the properties of the critical Casimir force in 
thin films of systems belongin g to the Ising universality class (UC) focusing on the crossovers 



between various surface UCs 



19 



20|, i.e., systems for which one or both surfaces give rise 



to relatively weak adsorption. Representatives of this class are simple fiuids, binary liquid 
mixtures, or Ising ferromagnets. The temperature dependence of the critical Casimir force, 
its sign, and its strength depend on the nature of the confining surfaces, which impose 
specific boundary conditions (BCs) on the relevant order parameter profile. In this context, 
so far only the cases of strongly adsorbing or neutral surfaces forming various surface UCs 
[s], is], [l4| have been studied theoretically, for both symmetric and antisymmetric BCs. In 
order to understand and thus to be able to control the aforementioned tunability of the 
critical Casimir forces, here we focus on the crossovers between these different surface UCs, 
i.e., systems for which one or both surfaces give rise to relatively weak adsorption. 

In order to calculate the critical Casimir force in the spirit of fieldtheoretical renormal- 
ization group theory we use the Landau- Ginzburg model in the film geometry. Within this 
approach the surfaces 1 and 2 are characterized by surface fields /iq,*; = 1)2, conjugated 
to the order parameter at the surface, and by so-called surface enhancement parameters 



co,i, i = 1,2, describing the tendency of the system to order at the surface [19 



Our 



results have been obtained numerically within mean field theory (MFT) as the lowest order 
contribution in a systematic 4 — d expansion in d spatial dimensions. The universal scaling 
functions of the order parameter profile and of the critical Casimir force have been calculated 
and thoroughly analyzed for the crossover between the so-called normal and the so-called 



special sur 
transition 



ace transition and for the crossover between the normal and the ordinary surface 



19 



20[, i.e., for various values of the parameters /io,i and co,j. It turns out that 



depending on the choice of these surface properties, the critical Casimir force can change its 
sign once or even twice upon varying the temperature. We also propose a simple criterion 
relating the sign of the critical Casimir force to the values of the order parameter at the 
surfaces and in the bulk. 

Recently, the crossover behavior in the same type of model, but for symmetry-preserving 
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BCs only (which enables one technically to go beyond MFT) has been studied by field- 
theoretic methods 2l|. Explicit two-loop renormalization group calculations show that the 
critical Casimir force can be of either sign depending on the surface enhancement parameters 
Co,i and co,2- However, these results are not applicable for fluid systems, because generically 
these are exposed to symmetry breaking surface flelds. Here we study the experimentally 
relevant case of tuning surface fields. In the presence of arbitrary surface fields both Monte 
Carlo (MC) simulation data 22| for the three-dimensional Ising model in the film geometry 



23 



25| show similar trends in the 



as well as exact results for two-dimensional Ising strips 
behavior of the critical Casimir force. 

Our presentation is organized as follows: In Sec. [IT] we introduce the model and briefly 
present the relevant basic theoretical facts concerning flnite-size scaling and surface UCs. In 
Sec. mil we report results of our calculations for the crossover between the special and the 
normal transition (Subsec. IIII A|) and the crossover between the ordinary and the normal 
transition (Subsec. IIII Bl) . In Sec. IIVI our results pertinent to four spatial dimensions are 
compared to results for d = 2 available in literature 23|-|25[ and we provide a comparison 
with MC simulation data for (i = 3 Ising fllms |22|. Section IVj summarizes our results. 



II. THEORETICAL BACKGROUND 



A. The model 



In order to calculate the order parameter proflle and the critical Casimir force we use the 
standard reduced Landau-Ginzburg Hamiltonian 1-1 = 1-1/ (ksT) (in units of the thermal 
energy fc^T and thus dimensionless) describing a system with 0{N) symmetry. For the fllm 
geometry with planar, laterally homogeneous surfaces and within MFT the order parameter 
(OP) proflle $ depends only on the spatial variable z orthogonal to the surfaces so that 
n = An with 

n (^)] = /'^^' I i ' + (z) + (z) - /.o,5<f (z) 

and where A is the macroscopically large, {d — l)-dimensional area of one of the equally 
sized conflning surfaces. Corresponding to the Ising UC studied here $ is a scalar. $i = 

4 



(1) 



^ {z = —L/2) and ^2 = ^ {z = L/2) are the values of the order parameter at the confining 



walls. The coefficient r oc 



where T^b is the bulk critical temperature, changes sign 



at bulk criticality. The coupling constant g > stabilizes for T < Tc^b- In the following 
we assume that the ordering bulk external field hbfl is zero, i.e., we focus on the critical 
concentration of the fiuid. The effects of the surfaces on the system are captured by the 
surface fields /io,j and by the surface enhancements Co,i as will be discussed in Subsec. IIIBi 
In the sense of renormalization group theory Eq. ([1]) captures all relevant scaling fields and 
thus is able to predict the leading universal behavior of critical films 
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2Q|. 



In the film geometry the critical Casimir force per area A of one of the equally sized 
confining surfaces and in units of ksT is given by 



c 



dL 



(2) 



where the excess free energy per area and in units of ksT is defined as 



r = {i-ib) L/{kBT). 



(3) 



Here f is the singular contribution to the total free energy of the film per volume V = LA 
and fb is the singular part of the bulk free energy density. 

Within MFT the analysis of the critical Casimir force leads to the Euler-Lagrange equa- 
tion (ELE) 



92$ 

dz'^ 



with the BCs 



(9$ 

dz 



r<l> + 

6 



co,i$i - ho^i 



(4) 



(5) 



-L/2 



and 



dz 



(Co,2$2 - /i0,2) • 



(6) 



z=L/2 



Instead of using the definition of the critical Casimir force given by Eq. ([2]), it is more 
convenient to use the thermal average of the {z, 2;)-component of the stress tensor 3|, |26| (we 
omit here the brackets (■) indicating the thermal average): 



- (z) , 



(7) 



which is spatially constant throughout the film including the surfaces. In order to obtain 
the critical Casimir force the bulk contribution has to be substracted: 



fc = [$] - , (8) 

where the bulk value is 

= ±e (-r) (9) 

and where G is the Heaviside step function. 

Dimensional analysis leads to the following scaled, dimensionless quantities: 

C = z/L (10) 

M (0 = yf ^^/^<^> (c^) {p/'^)mft = 1 (11) 

y = rL^'" vuFT = 1/2 (12) 

Ci = Co,iL<^/'^ (0/^^)mft = 1 (13) 

= ^Jlho.L^'^"/^ {KV^)mft = 2, (14) 

where is the crossover exponent and A^^ is the surface counterpart of the bulk gap exponent 
A. In terms of these scaled variables the Hamiltonian in Eq. ([1]) takes the scaled form 

n (.)] J-L-^ r I i + ',M^ (0 + W (0 



9 J -1/2 ^ {2\ dC J 2 ^^^4 J (^5) 

+ y-^? - hiMi + ^Ml - /12-M2, 

where within MFT the prefactor 6/(7 is undetermined. In these scaled units the bulk limit 
(which minimizes the integral in Eq. ( IT5|) for constant M. = Aih) and the correlation length, 
defined via the exponential decay of the two-point correlation function for T 7^ Tc,6, are 

Mb = ±Q{-y)\y\''' (16) 

and 

respectively. 
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The above MFT considerations can be put into the general context of scahng theory. 
In accordance with this the order parameter profile in a film exhibits the following scaling 
behavior [l^: 

<!>{z,t,L) = 'i>,{t)M{z/^,L/0, (18) 

where ^ (t = (T — Tc^h) /Tc,b — ^ ±0) = I'^^r'^; are nonuniversal amplitudes with a uni- 
versal ratio {it / prp = "\/2- ^b{t) = '^'b,o l^l'^ is the bulk order parameter (we take 
$6 {t > 0) to be $6 (— t) valid for T < Tc^t) with the (only) second nonuniversal amplitude 
^b,o- In Eq- ( fTSj) Ai {x, y) is a universal scaling function which is normalized such that 
M. [x oo,y = oo) = 6 (— t). The scaling function Jvl {x,y) is suitable for capturing the 
semi-infinite limit. For finite film thickness it is helpful to rewrite Eq. f lTS]) in terms of a 
different, also universal scaling function Ai = z/L,y = (sign t) {L/^Y^" = t [^L/^qY^^^: 

$ (z, t, L) = % (t) i^/Lf'^ M [z/L, t (L/U) '/^) = $,,0 (L/^^) '^'"^ M (C, y) . (19) 

Note that in Eq. f lT9|) we have used as expression for ^ the one above also for t < 
so that here ^(t) = ^o'|t|~'^ for all t; this yields a scaling variable y which is ana- 
lytic in t. The universal scaling function Ai {(, y) has the property and is normalized 
such that, for Clz/T = const, (C — )• 0, y — )■ oo) = Q {—y) \yf = M.h{y) which implies 
lim \{i/Lf''' M (z/L, (sign t) {L/^)^^A] = 6 {-t). In the following Mb{y) will be also 

t fixed 

called as bulk order parameter. The above MFT expressions (Eqs. (lT0l)-([l4l)) are in line 
with these general properties by noting that within MFT one has r = (.^J") ^ t. 



Finite-size scaling theory [19|, |20|, |27|, |28| and renormalization group theory 



29 for the 



film geometry predict that the critical Casimir force takes on the following scaling form: 

fc {-T, Co,i, /io,i, Co,2, ho,2; L) = {d- 1) (y, ci, hi, ca, /la) • (20) 

Within MFT the universal scaling function d of the critical Casimir force can be determined 
only up to an undetermined prefactor ~ 1/(7 (compare the note after Eq. ( fTSl) ). An appro- 
priate way to cope with this is to express d in units of the critical Casimir amplitude A(+^+) 
at Tc^f, for fixed point BCs hi = h2 = oo, which carries the same undetermined prefactor. 
Accordingly, all our MFT results expressed in units of A(+^+) {d = 4) = — {6/g) < 0, 



where K = K{l/2) is the complete elliptic integral of the first kind 30|, are independent 



of this undetermined prefactor and are therefore accessible to comparisons with results ob- 
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tained by other theoretical techniques or experimentally. (Note that within MFT, i.e., d = 4, 
g is dimensionless.) 



B. Surface universality classes 

Near a surface the system properties differ from their bulk values. For example, the 
tendency to order is reduced due to missing neighbors and the coupling between the ordering 
degrees of freedom at the surface can differ from its bulk value. In Eq. ([T]) these surface 
effects are captured and characterized by the surface enhancements Co,i, which preserve the 
0{N) symmetry. In addition, the surface can favor one bulk phase over the other. In Eq. ([1]) 
this explicit breaking of the symmetry is described by the surface fields /io,i- For detailed 



discussions of surface criticality see Refs. 



19|, 



2Q|, 



31| and Ref. [321]. In the latter, microscopic 



expressions for ho and Cq are derived by using density functional theory. 

A semi-infinite system the surface of which has a reduced tendency to order belongs to 
the so-called ordinary surface UC (in the following labeled by u = ord). The corresponding 
fixed point values of the surface parameters are Cq^'^ = oo and h^^'^ = 0. Because cq is a 



20|, 



so-called dangerous irrelevant variable, it may not simply be set to its fixed point value 
but rather the following linear scaling field has to be considered: 

hr = ho/cl (21) 

with the scaling exponent 

=1. (22) 



Y MFT 

Thus the dependences on the scaling variables in Eq. (EU]) reduce to the dependence on a 
single scaling variable per surface, which, within MFT, is V^f^ = hi/ci ~ {hQ^i/co^i) L^°"^/^ 
where Af^ = 1 /2 within MFT. 

Systems with surfaces, which do not break explicitly the 0{N) symmetry of the order 
parameter but locally enhance the ordering, belong to the so-called extraordinary surface 
UC characterized by the fixed point values Cq^ = — oo and h^^ = 0. Such surfaces order at a 
surface transition temperature Tc^s > T^^b- In the parameter space spanned by (T, cq, = 0) 
the transition lines of the ordinary transition, the extraordinary transition, and the surface 
transition meet at the multicritical point {T^^b, Cq) of the so-called special transition (in the 
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following labeled by u = sp). Within MFT the value of the surface enhancement at the 
special transition is Cq = 0. Near the special transition both ho and cq are relevant. 

Surfaces which prefer one of the bulk phases over the other break the symmetry of the 
order parameter explicitly. Semi-infinite systems with such surfaces belong to the so-called 
normal UC (in the following labeled by u = ±, depending on the sign of the surface field) 
characterized by the fixed point /iq°'~™ = oo. For both normal and extraordinary UCs the 
surface is ordered even above Tc^b- Explicit calculations for the Ising model |33f] and renor- 
malization group analyses |3l|| reveal the equivalence of the normal and the extraordinary 
transitions. 

For films the character of both surfaces matter. Therefore the UCs for the systems 
in the film geometry can be labeled by the two UCs of the corresponding surfaces in the 
semi-infinite geometry. This interpretation is meant implicitly when in the following the 
crossover behavior in the film geometry is named after the crossover between surface UCs 
in the corresponding semi-infinite systems. 

At the fixed points and within MFT the scaling functions 'i9ui,u2) ^ sp,ord}, are 



known analytically [30|, |3^. In order to obtain the scaling functions for the surface param- 
eters Co,j and h^ i being off their fixed point values, we have minimized the Hamiltonian in 
Eq. ([1]) numerically. 



III. RESULTS 



A. Crossover from the special to the normal transition 

In this subsection we study the behavior of the scaling functions for various applied surface 
fields /io,i and /io,2, including /lo.i 7^ ho,2, at vanishing surface enhancement parameters 

Co,l = Co,2 = 0. 



1. Rescaling 

For a wide range of finite values of the surface fields hi, the variation of the scaling 
function {y) of the critical Casimir force as a function of the scaled temperature y (see 
Eq. ( !T2|) ) is very similar to the one corresponding to the fixed point solutions ^{+,+) (y) or 
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i.y)i depending on the sign of hih2 (see Fig. [T]). As described below this similarity can 
be specified quantitatively by a suitable rescaling of the scaling functions. 

This rescaling idea is borne out by our observation that the OP profile $ {z\ r, /io,i, /io,2, L) 
of the film of width L exposed to finite surface fields /lo.j, "^ = 1,2, can be expressed in terms of 
the OP profile $ [z] r, /io,i = oo, /io,2 = ±00, L*) = $(+,±) {z; r, L*) of a film of width L* > L 
with infinite surface fields, in such a way that the former profile is a portion of the latter 
one. (We recall that for films with infinite surface fields the OP diverges at the surfaces.) 
Since for both scaling functions of the OP profile the spatial scaling variable has the range 
—0.5 < C ^ 0-5, the scaling function corresponding to the film with finite surface fields has to 
be shifted and "stretched" in order to fit into the corresponding interval. The corresponding 
scaling variable and the amplitude of the scaling function may be different from those of 
•^(+,±) (C; y)- Therefore we make the ansatz 

M (C; y, h, h^) = fM(+,±) (r (c - Co) ; y*) • (23) 

Comparing the ELE for both scaling functions, i.e., the scaled form of Eg. (H), Ad" = 
yM+M^ and A^('+,±) = l/*-M(+,±) + -^(+,±). where M' = dM/dC, we find |35| 

/ = r, y* = r-^y. (24) 

The rescaling function r = r (y, /ii, /12) and the shift Co = Co (y, ^1)^2) are obtained 
from the scaled form of the BCs in Eqs. (jSD and dHD, —hi = A^' (— 0.5; /ii, /i2) and 
/12 = -M.' (0.5; y, hi, /i2), respectively: 

^'-^(+,±) {r (-0.5 - Co) ; r-^y) = -hi (25) 
^'-^'(+,±) (0-5 - Co) ; r-^y) = h2. (26) 

We note that Co 7^ only if \hi\ ^ |/i2|. With this the corresponding mapping for the 
scaling function of the critical Casimir force can be read off from the expression for fc in 
terms of the stress-tensor, i.e., from Eqs. ([7]) and ([H]) by using Eqs. ffTT]) . flTB]) . and fl2U]) 
(note that Eqs. ([7]) and ([8]) are valid only within MET and thus require ci = 4 so that g is 
dimensionless) : 

^ {y- hi, h,) = {d- 1)-' ^ Q (Mf -1{M'- M,') -1{M'- M,')^ , (27) 
from wich we find 

^iy;hi,h2)=r'^i+,±){r''y). (28) 
10 



The ranges of values of the surface fields hi, for which the proposed mapping can be ap- 
plied, is limited by the necessity to fulfill Eqs. ([25]) and (ES]). For the (+, +) BC A^(+,+) (C; v) 
varies from — oo to for —0.5 < C ^ and from to oo for < ( < 0.5. Thus, for any 
hi > and h2 > Eqs. f l2^ and f l2B]) have a solution for r and Co- We have checked the pro- 
posed mappings (Eqs. f l23p and (128 p ) by comparing them with the results of the numerical 
minimization of the Hamiltonian in Eq. ([1]) and we found excellent agreement. Examples 
for the mapping of are shown in Fig. [2j 

In the case of opposing surface fields, —hi<h2<0 (here \hi\ > |/i2|, otherwise 1 and 
2 have to be interchanged), the mapping is restricted by the upper bound for the slope 
■M^ I (C; y) < 0. For fixed (but arbitrary) hi > 0, only films with —hi < h2 < h2^min < 



361] can be mapped onto films with (+, — ) fixed point BC, where h2^min depends on both hi 
and y. For y > 0, h2,min {hi, y) is obtained by evaluating A^'(_|_ _) in Eq. fl26|) at the position 



where 



is minimal, i.e., at Co = 0.5 and we find 

^2,mm {hi,y) = rl^inM\^^_) (O; r^^„y) (29) 
with Tmin given implicitly by Eq. ( 125|1 : 

^mm-^(+,-) (-'"mm;r--„l/) = ~hi. (30) 

Since A^(+,_) (C) is defined for -0.5 < C < 0.5, Eq. ^ together with Eq. ^ yield the 
condition < r^m < 0.5; r^m attains its limiting value 0.5 for hi = oo independent of y > 0. 
\h2,min\ is largest for ?/ = and decreases upon increasing y. For y < the determination of 
h2,Tnin is more involved, because A^(_,_ „) is no longer monotonic for —0.5 < C ^ 0. For fixed 
hi, the maximum of {} {y; hi > 0, /i2 < h2^min {hi,y = 0)) as a function of y moves towards 
y = upon decreasing /z.2 and we observe that for /i2 = h2^min{hi,y = 0) the maximum 
occurs at Tc^b- 

It turns out, that for sufficiently strong surface fields the rescaling functions r {y, hi, h2) 
and Co {Vi hi, h2) are well approximated by resorting to the short distance behavior of M. and 
they become rescaling parameters which are independent of y. One obtains (see Appendix |A|) 

T,d=[l + l{hi) + l{h2)]-' (31) 

and 

Q'={l{h2)-l{hi))/2, (32) 
11 



where I [hi) = 2^/^\hi\-^^\ The comparison with our numerical data shows, that for systems 
with hi,h2 > 10 this approximation works very well for the whole range of y: for these films 
the difference between using r {y; hi, /12) and Vsd {hi, /i2) for the mapping is not visible on the 
scale of Fig. [2l On the other hand, for values hi, —h2 > 10 this approximation is only valid 
for y > ymax, where ymax is the position of the maximum of {}. As expected, our numerical 
data show, that at T^^b for those values of the surfaces fields for which the short-distance 
approximation holds, the OP profile of the film near the surfaces can be well aproximated 
by the corresponding OP profiles of the semi-infinite sytem. For weaker surface fields the 
profiles differ significantly - even at the surfaces. 

The deviation of the scaling function {y; hi, h2) of the critical Casimir force from the 
fixed point scaling function ^(+,±) {y) in leading order of hi and /i2 is obtained from the 
mapping given by Eq. (1251) and from r^^ in Eq. (IHT]) : 

^{y;hi,h2) -%,±) (y) ^ 

- 2^/4 {\hi\-'^' + \h2\-'^') (y) - (y)] , \hi\ , \h2\ ^ 00. (33) 

Our numerical data for the scaling functions are in agreement with Eq. ( 133|) for , |/i2| ^ 
1000. This algebraic behavior explains the slow convergence of the scaling function i} towards 
the fixed point scaling function '&(+,±) as apparent from Fig. [TJ We note that, since the term 
in curly brackets in Eq. (!33|) is comparable with ^(+,±) itself, the relative deviation for, e.g., 
hi = 1000 and /;,2 = C)0 is still about 15%. Our data for the critical Casimir amplitude 
A+ (hi) = A {hi, /i2 = hi) = {} (y = 0; hi, /i2 = hi) shown in Fig. |3]^a) display an algebraic 
behavior which is in accordance with Eq. fj33|) . 



2. Weak surface fields 



For weak surface fields /io,i = /io,2, i-e., if the length scale Isp = (^\/ 9 /QhQ i^ 37| 
associated with the surface field in the semi-infinite geometry dominates over L, both surfaces 
approach the special transition (we recall that here we consider the case Ci = C2 = 0) and 
the order parameter profiles do not vary substantially across the film (data not shown). In 
this case the square gradient term in (Eq. ([T])) can be neglected and the order parameter 
can be approximated to be constant so that the unsealed free energy (per fc^T and A) is 

f iho,i = ho,2 ^0) = L (^$2^^ + |$to..) - 2 ho,i (34) 

12 



Comparison with Eq. ([T]) shows, that the surface fields /io,i = ho,2 act hke an effective 



bulk field hb^Q = 2ho^i/L. This is in line with the analysis by Nakanishi and Fisher [38| 
of the shift of the critical point and of the phase boundary in films, which shows that for 
^0,1 = ^0,2 > the phase boundary is shifted towards negative values of the bulk ordering 
field /io,fe- Minimization of the free energy in Eq. with respect to i yields at r = 
the critical value of the order parameter: 

/ 12 X 1/3 



= ■ '''' 

The dependence on /io,i, i.e., the exponent 1/6 = 1/3, is the same as the one with which 
at Tc^b the bulk order parameter responds to a weak bulk field; Eq. ( l35i) agrees with our 
numerical data (not shown). Accordingly, the free energy (per kf,T and A) at criticality is 

fcr {ho,l = ho,2 ^ 0) = i^T^) ■ ^^^^ 

Therefore the critical Casimir amplitude A+ (hi) depends on the surface fields as 

A+ ih) ~ -ht^\ hp > L. (37) 

This checks with our numerical data shown in Fig. |3]^b). Note that in the regime L < Igp the 
film thickness has not yet reached asymptotically large values. As a consequence fc ~ L~^/^ 
which is a much slower decay than ~ L~'^ in the asymptotic regime L ^ l^p (Eq. ( |20|) ). 

The change from a purely repulsive scaling function of the critical Casimir force with a 
maximum below T^.^ (which are the characteristics of the (+, — ) BCs) to a purely attractive 
scaling function of the critical Casimir force with a minimum above Tc^b (which are the 
characteristics of the (+, +) BCs) occurs at fixed /ii > for h2 — > 0~, i.e., in the regime 
where the rescaling scheme is not applicable. We note that in this crossover regime the 
overall magnitude of d is much smaller than the one for the (+, — ) and (+, +) BCs (for all 
curves shown in Fig. [T]the rescaling scheme does apply). As an exemplary case, we show 
in Fig. m d {y) for hi = 168 and various values of /i2- For ^2 = —4.2 ^ h2^min (2/ = 0), so 
that h2/hi = —0.025, the maximum of i} is located at ?/ ~ 0. Upon a further decrease of 
|/i2| a minimum develops at some positive value ymin and exhibits two broad maxima of 
comparable height, one above and one below Tc^b- We have found, that the value of the 
scaling function at the maxima as a function of /i2 varies as (ymax) ~ (^2)^ for ^2 — ^ 0. For 
a further decrease of |/i2| the minimum becomes negative. For all /i2 < the scaling function 
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present two different approaches to describe analytically the variation of {y; hi, in the 
crossover regime. These approaches deal with films for which the rescaling scheme is not 
applicably and within suitable ranges of surface fields they capture well the numerical data 
for the crossover behavior shown in Fig. |H 

B. Crossover from the ordinary to the normal transition 

A surface with a nonzero surface enhancement co,i > and with /iq,! 7^ falls into 
the crossover regime between the ordinary and the normal transition. From the boundary 
condition in Eq. ([5]) it follows that within MFT for a surface field /io,i = Co,i$(, the semi- 
infinite profile is spatially constant: $ {z) = Here we consider only the case /io,6 = 0^ 
(see Eq. ([T])) for which ± coincides with the sign of /io,i (i-e., the phase preferred by the 
surface is the same as the one prevailing in the bulk) in order to avoid complications induced 
by wetting transitions which occur if the surface preference is opposite. Thus for stronger 
(weaker) surface fields the value of the order parameter at the surface is larger (smaller) 
than in the bulk. For r > 0, = and the surface is ordered for any h^ i > 0, resembling 
the normal transition. On the other hand, from Eq. it follows that for temperatures 
sufficiently below Tc^b, i-e., 

r < r^,i = - (5-/6) (/io,i/co,i)% (38) 

where r^^i is given implicitly by (r^) = /io.i/co,i (which happens to be the wetting tran- 
sition temperature at wall 1 for /io,6 = 0^ [401]), due to the suppressive influence of co,i 
the order at the surface is lower than in the bulk and increases with increasing distance 
from the surface, resembling the order parameter behavior corresponding to the ordinary 
transition. Therefore, for fllms exposed to nonzero and finite surface parameters we can 
identify the following regimes resembling different surface UCs: (1) far above Tc^f, the (+, ±) 
UC (depending on the relative sign of the surface fields) and (2) far below T^^ (Eq. (1551) ) 
the {ord, ord) UC. For temperatures in between, a regime resembling the (+, ord) UC can 
exist. The crossover between these different regimes gives rise to richly structured scaling 
functions of the critical Casimir force (see below). 
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Since the critical Casimir force is repulsive for films with (+, — ) and (+, ord) BCs and 
attractive for films with (+, +) and {ord, ord) BCs, in view of the above discussion the critical 
Casimir force is expected to change its sign up to two times. Let us first consider the case 
^o,i/co,i > ^o,2/co,2 > so that t^^i < t^^2- Accordingly the change from an attractive to a 
repulsive force is expected to occur at a certain temperature Tc^2 = Tw,2 + ^T2,i upon lowering 
the temperature (crossover from the regime resembling the (+, +) UC to the one resembling 
the (+, ord) UC), followed upon further lowering the temperature by a change from repulsion 
to attraction again at Tc^i = r^^i + Sti^2 (crossover from the regime of the (+, ord) UC to the 
one of the {ord, ord) UC). If /io,i/co,i > — ^o,2/co,2 > only one change from a repulsive to 
an attractive force is expected to occur at the temperature r^i = r^^i + 5ri 2 (crossover from 
the regime resembling the (+, ord) UC to the one resembling the {ord, ord) UC), because at 
Tc,2 = T«,,2 + 5^2,1 the film crosses over from the regime of the (+, — ) UC (due to /io,i^o,2 < 0) 
to the one of (+, ord) UC with both UCs rendering the critical Casimir force to be repulsive. 
6Tij are the finite-size corrections depending on the properties of both surfaces. The changes 
of sign and the exact values of r^i will be discussed in more detail in Subsec. IIII CI 

In the following we restrict ourselves to a more qualitative description of the features 
in the crossover regime, although for particular regimes the rescaling and the perturbation 
theory, as introduced in Subsec. UlI A H and in Appendix[Bl respectively, should be applicable 
as well. 

First we discuss the symmetric case ci = C2 and hi = h2. In Fig. [5] we show our 
numerical data for the scaling functions of the critical Casimir force "dsymm {y, 100, hi) = 
^ {y; ci = C2 = 100, hi = /i2). We note that ^symm is negative for all values of y as in the 
case of the (+, +) and {ord, ord) BCs. It exhibits a minimum above Tc^, like the fixed point 
scaling function ^{+,+) [30]. This minimum is very shallow for weak hi and, as expected, 
deepens for stronger surface fields hi. As discused above, for a more negative scaling variable 
y the film crosses over to the asymptotic regime of the {ord, ord) BCs. This results in the 
appearance of an additional minimum below Tc^, as it occurs for the scaling function ^[ord,ord) 
34| . This minimum deepens for decreasing hi. For finite and nonzero values of the surface 



parameters {ci,hi) the cusp-like minimum of ^(ord,ord) (which is a MFT artefact 34j. l4l|) 
is smeared out (see Fig. [5]). For ci = 100 the two minima are equally deep for hi = 168. 
Between the two minima ^symm exhibits a maximum at the value y = yw,i, corresponding 
to Tw,i (Eq. (l38l)). with ^symm {yw,i) = for all hi = h2 and Ci = C2 > 0. (For very weak 
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(very strong) hi the minimum above (below) Tc^ is very shallow and finally disappears for 
/ii — 7- {hi — 7- oo) and thus also the maximum in between becomes hardly detectable.) 

For general values of the surface parameters, the crossover values ?/c,i and ?/c,2; correspond- 
ing to Tc^i and Tc^2, respectively, differ and, as discussed above, the film crosses through an 
additional regime corresponding to the (+, ord) BCs. This is illustrated in Fig. E] for the 
case ci = C2 = 100, hi = 168, and different /i2- 

For large surface fields with opposing sign, hih2 < 0, the system is in the asymptotic 
regime of the (+, — ) BCs, i.e., the critical Casimir force is repulsive and exhibits a maximum 
below Tc^b (see the curve for h2/hi = —4 in Fig. [H]). With decreasing \h2\, according to 
Eq. fl38p the asymptotic regime of the ordinary transition is reached at the second surface 
already for less negative values of y: the scaling function of the critical Casimir force becomes 
negative at negative values of y and exhibits a minimum there (see the curves for h2/hi G 
{—2, —0.9} in Fig. [6]). For the case of fixed Ci = C2 = 100 and hi = 168 the minimum is 
deepest for h2 — —0.9hi. 

For /i2 = the scaling function of the critical Casimir force is positive for y > ?/c,i — —3.6 
and exhibits a maximum (ymax) / |^(+,+)| — 0.025 at ?/c,i < ymax — — 1 < 0, like for the 
(+, ord) BCs (see Fig. [6]). For y < the film is in the asymptotic regime of the {ord, ord) 
BCs and the scaling function of the critical Casimir force is negative with a minimum 

^ {ymin ^ -14) / |A(+,+)| ~ -0.114. 

For hih2 > and for large positive and very negative values of y the film is in the 
asymptotic regime of the (+, +) BC and the {ord, ord) BC, respectively, and thus is 
negative in these limits. However, for large values of h2, \yc,2\ becomes very large such that 
the limit of negative is not seen anymore (see the curve for h2 = lO/ii in Fig. |6]). In the 
intermediate regime between ?/c,i and ?/c,2, corresponding to the (+, ord) BCs, ^ is positive 
with a maximum at negative y (see the curves for h2/hi G {2, 10} in Fig. [6]). 

Close to the ordinary transition there is only one relevant, linear scaling field h"'"'^ = hi/ci 
associated with a single surface (Eq. (I2T]) ). As expected, we see in our data that the scaling 
functions -(9 {y) for different values of q and hi but fixed h°^'^ = hi/ci indeed exhibit similar 
behaviors. For small Cj and hi, the details of the shape and the amplitude of still depend 
on their particular values, but with increasing values of these two surface parameters a 
convergence of the corresponding i) is seen. 

In order to illustrate this point, in Fig. [7] we discuss the case in which C2 is varied at fixed 
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Ci = 100 and hi = h2 = 168. For C2 > and large positive (negative) y the film corresponds 
to the (+, +) BCs (the {ord, ord) BCs) and is negative with a minimum above (below) 
Tc^b (for C2/C1 > 50 partly not visible on the scale of the figure). In between there is an 
intermediate regime, corresponding to the (+, ord) BCs where is positive with a maximum 
below Tc^fj- With increasing C2, i.e., surface 2 approaching the ordinary transition, becomes 
positive already at larger values of y and finally the regime of negative for positive y 
becomes hardly visible (see the curve for C2 = 50ci). On the other hand, the minimum at 
negative y deepens with increasing C2. We point out, that this is the same qualitatively 
behavior as for fixed Cj and /i2/^i 0^ (compare the discussion above and Fig. [H]). For 
decreasing C2 the minimum at positive y, corresponding to the (+, +) BCs regime, depeens 
(compare the curves for C2/C1 G {5,1,0.5,0} in Fig. [7]). The regime corresponding to the 
{ord, ord) BCs with negative finally disappears for C2 — )■ (see Fig. [7]). The same behavior 
is seen for increasing h2, since it leads to a stronger ordered surface as decreasing C2 does 
(compare the curve for h2/hi = 10 in Fig. [6] and the curve for C2/C1 = in Fig. [7]). 

For a negative surface enhancement C2 < the trends of deepening of the minimum above 
Tc^b and increasing of the maximum below Tc^ continues (see the curve for C2/C1 = —0.1 in 
Fig. E]), because the extraordinary transition, which corresponds to C2 < 0, and the normal 
transition are equivalent (see Subsec. Ill BO . 



C. Change of sign of the critical Casimir force 

In the crossover regime we observe the interesting feature that the scaling function of the 
critical Casimir force changes sign as function of the scaling variable y = t {L/^qY^'^ ■ This 
implies that at fixed temperature the critical Casimir force changes from being attractive to 
being repulsive (or reverse) upon varying the distance between the two surfaces. Equivalently 
this change also occurs at a fixed distance upon varying the temperature. Thus temperature 
allows one to control both the strength and the sign of the critical Casimir force. 

In order to understand and to interpret this change of sign we consider the functional 
form of the order parameter profile ^semi {z > ~Lj1\ r, ©1) in the semi-infinite system with 
the surface located at 2; = —L/2 and with the surface parameters ©1 = (co,i, /io.i); this is 



available in the literature, e.g., in Ref. 



42|. Using the BC (for z = L/2) in Eq. for 



certain values of the parameters ©2 = (co,2, ^0,2) of a second surface located at z = L/2 one 
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can find a temperature t = Tc such that 



ho,2 = Co,2$.emi (z = L/2; Te, ©i) + ^'^^^^ (z = L/2; r„ ©i) . 



(39) 



Upon construction, at this temperature Tc, ^semi [z] t^, &i) coincides with the profile in a 
film (of width L and with confining walls ©i and ©2) for distances —L/2 < z < L/2. 
Equation fl39l) implicitly defines classes C {L,Tc] &i) of walls ©2 which do not disturb the 
semi- infinite profile if they are inserted at the distance L and at the temperature r^. Since 
such OP profiles of finite and semi-infinite systems coincide, the critical Casimir force for 
a film of thickness L with surfaces ©1 and ©2 G C{L,Tc]&i) at this temperature Tc is 
zero. Because in the semi-infinite geometry the order parameter varies algebraically for 
\z + L/2\ and decays exponentially towards its bulk value for \z + L/2\ ^ ^ (see, e.g., 



Ref. [3l|), for sufficiently large L the value of the order parameter at the position of the 



second "nondisturbing" surface ©2 G C (L, ©1) and for the temperature Tc is close to its 
bulk value. Thus, in a first approximation, the change of sign occurs if at one of the surfaces 
the order parameter takes its bulk value, = We recall, that a surface which enhances 
the order, i.e., at which |$j| > corresponds to the normal transition (±) while a surface 
that suppresses the order, i.e., at which |$j| < \^b\, corresponds to the ordinary transition 
(ord) and that the critical Casimir force is attractive (repulsive) for (+, +) and {ord, ord) 
BCs (for (+,-) and (+,ord) BCs). 

Based on these features, within MFT the sign of the scaling function -(9 of 
the critical Casimir force can be inferred rather reliably from the values Aii = 
M {( = —0.5; y, Ci, hi, C2, /i2) and A^2 = (C = 0.5; y, Ci, hi, C2, /i2) of the scaled order pa- 
rameter at the surface 1 and 2, respectively: 



In Fig. [H] this is illustrated for two examples. The case Ci = C2 = 0, hi = 2.5, and 
/i2 = — O.l/ii is representative for the crossover regime between the normal and the special 
transition, like the cases shown in Fig. IH note, however, that there hi = 168 so that 'd differs 
even if the ratio h2/hi is the same. The scaling function ^ {y; Ci = C2 = 0, hi = — 10/^2 = 2.5) 





(40) 
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is positive for |?/| ^ 5 and decays to zero for |?/| — )■ oo (see Fig. [8]^a)). It exhibits maxima 
at ymax,i - -4 and ymax,2 ^ 9 of approximately equal height -d {ymax,i) / |^(+,+)| - 0.00066. 
Between these two maxima d decreases towards negative values and exhibits a minimum 
•§ {ymin — 1) / |^(+,+)| — —0.0049. As shown in the inset of Fig. [H](a), for ?/ > 9 the relation 
A^i > = > 0) > M.2 holds and the critical Casimir force is repulsive. For ?/ ~ 8, /i2 

is too weak to overturn the positive order imposed by hi and M.2 turns positive. For y < S 
the criterion proposed in Eq. fHOl) predicts an attractive critical Casimir force, whereas the 
actual change of sign occurs at y ~ 5. This discrepancy occurs because the critical Casimir 
force is zero if $2 = ^semi {z = L/2, r, ©1) whereas in Eq. fHU]) . which is proposed to hold 
for large L, we assume ^semi {z = L/2, r, ©1) ~ Below Tc^b, the bulk value of the order 
parameter is nonzero (and positive, see Eq. flTB]) and the horizontal dotted lines in the insets 
of Fig. [8]). For decreasing, i.e., more negative values of y, the increase of the positive order at 
the second surface (c2, /i2) is reduced by the opposing surface field /i2 < 0, whereas at the first 
surface (ci, hi) the positive order is enhanced. At approximately y —3 the bulk order Aij, 
and the surface order A^2 become equal and for y ^ —3 the relation 0<A^2<A^b<A^i 
for the order parameter holds. Hence, according to Eq. ( HOj) a repulsive critical Casimir 
force is expected. As can be seen in Fig. M^a) this is indeed the case, but the change of sign 
occurs at y ~ —1.5. Again, this discrepancy occurs due to the asymptotic character of our 
simplified criterion in Eq. PUI) . 

As a second example Fig. El^b) shows the scaling function of the critical Casimir force 
for the case ci = C2 = 100 and hi = 10/i2 = 1680 corresponding to the crossover between 
the normal and the ordinary transition. The scaling function of the critical Casimir force is 
positive for y < —1.5 and negative otherwise. In general, for strong surface parameters c and 
h considered here the values of the order parameter at the surfaces do not change significantly 
upon varying the scaling variable y (compare the inset of Fig. [8](b)). For 50 > y > —50 
the values of the order parameter at the surfaces are in the range 14.9 < Aii < 15.5 and 
1.5 < A^2 < 2. As in the first example shown in Fig. [H](a), the change of sign of the scaling 
function of the critical Casimir force occurs approximately at that (here single) value of 
y for which Aih{y) equals A^2 (see Fig. [8](b)). Due to the strong surface field hi, Aii is 
larger than Ai^ for the whole range of y within which the scaling function of the critical 
Casimir force has a detectable amplitude. Far away from criticality the correlations reduce 
to microscopic length scales and the critical Casimir force becomes vanishingly small. Thus 
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in the case discussed here no attractive critical Casimir force is observed for very negative 
values of y. 



IV. COMPARISON BETWEEN d = 2, d = 3, AND d = A 



d = 2 and d = 4 



Recently 23|-|25| the crossover behavior of the critical Casimir force has been studied for 
a two-dimensional Ising strip in the limit of M ^ 1, where M is the number of rows of 
the strip, i.e., its thickness. These authors considered the case of an unchanged coupling 
constant Jii in the two surface rows of the strip. 



J|I = J(1 + A||), A||=0, 



(41) 



where J is the nearest-neighbor coupling constant in the bulk. Carrying out a systematic 



continuum limit one can relate the lattice parameters and t 



19 



model. The surface enhancement cq is related to A 

co= (l-2(d-l)A||)/ 
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le couplings in the continuum 



42|: 



(42) 



where a is the spacing of the simple cubic lattice. For Ay < (2 (c? — 1))~^ (within MFT) 
the continuum limit a — )• leads to Cq = oo which corresponds to the ordinary transition. 
The scaling variables x = sign{t) \t\'^ M/^q and z ~ h^'M, with u{d = 2) = 1 and 
A"^'' {d = 2) = Ai{d = 2) = 1/2, as introduced in Ref. 25| translate into the the scaling 
variables we use here according to x = y'^ and z = (h^'^'^Vj^^^ (compare Eqs. f|T2|) and f|T9|) 
and the expression after Eq. fl22l) ). Figure 9(a) in Ref. 25(] corresponds to the symmetric 
case and thus corresponds to Fig. [5] here. The qualitative behavior is the same in both 
cases. The scaling function i) (called y in Ref. 25|) is negative for the whole range of the 
scaling variable y. For weak surface fields, exhibits a minimum below Tc^f, which becomes 
more shallow upon increasing hi, while a minimum above Tc^b develops concomitantly and 
the absolute value of its depth increases. In contrast to d = 4, where the amplitude of the 
maximum between the two minima vanishes {i) [ymax] d = A) = 0), in d = 2 the amplitude 
of the maximum is nonzero (see, e.g., the dash-dotted curve for z = 0.5 in Fig. 9(a) in 



Ref. |2^). Concerning the antisymmetric case in d = 2 (see Fig. 9(b) in Ref. [25l]) and d = A 
(available data for /ii > and /i2/^i = —1 are not shown, but the qualitative behavior can 
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be inferred from the curves for /i2//ii G {—0.9,-2} in Fig. [6]) has one maximum with 
'^max > and one minimum with i)min < 0. In both cases, the minimum shifts towards more 
negative values of y for increasing hi. For d = 2 the position of the maximum is located at 
positive y for small hi, moves to smaller values of y with increasing hi, and for a particular 
value of hi it is located at y = 0. In contrast, for d = 4 and within the ranges of values of Ci 
and hi for which the value of the maximum is larger than numerical errors, the maximum is 
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25| the authors used also different scaling 



always located at negative values of y. In Refs. 
variables. We have compared the resulting scaling functions expressed in those variables, 
too, and have also found the same trends in d = 2 and d = 4. 

In Ref. 2J] and in its extended version 25| only the cases /12 = ±/ii were considered. 



Results for the critical Casimir force in the two-dimensional Ising film subject to arbitrary 
surface fields and for Ay = (Eq. fHTl) ) are provided in Ref. 3]. All three studies for d = 2 
show, that the critical Casimir force can change sign not only by changing the surface fields 
but also by varying the temperature. In the crossover regime the scaling function of the 
critical Casimir force exhibits more than one extremum and thus the behavior in ci = 2 is 
in qualitative agreement with the present MFT results. 



B. d = 3 and d = 4 



Due to a dearth of analytical means, the natural choice for studying the critical Casimir 
force in d = 3 consists of Monte Carlo (MC) simulations for a three-dimensional Ising film 
on a simple cubic lattice with (100) surfaces 43|]. Within this approach, the surface fields are 
mimicked by additional layers in which all spins are set equal to 1 or —1 and which couple 
to the actual surface layers of the film by a modified coupling constant Jj_ = a J, with J 
as the nearest neighbor coupling constant in the bulk {a = J±/J measures the modified 
coupling constant J± in units of J, not to be confused with the critical exponent a). For 
details concerning the determination of the critical Casimir force via such MC simulations 



see Refs. 
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43(1 . The universal critical Casimir amplitude for the d = 3 Ising film is 



A(_(_ (d = 3) ^ 0.38 [44]. The universal critical exponent of the correlation length for 
the d = 3 Ising bulk UC is u {d = 3) ^ 0.63 [45] and its non-universal amplitude for the 
simple cubic Ising model with nearest neighbor coupling is ~ 0.50 [46]. First results 
for a three-dimensional Ising film exposed at one surface to a fixed, finite surface field hi 
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corresponding to ai = 1 and at the other one to a variable surface field /12 corresponding to 
< 0:2 < 1 are shown in Fig. [91 These results show a similar behavior as the one discussed 
for d = 4 in Subsec. IIIIBt For h2/hi = 1, i.e., 0:2 = 1 the scaling function exhibits only one 
minimum which is located above T^^b- With decreasing h2/hi (i.e., smaller 0:2) this minimum 
becomes shallower and below T^^b a maximum develops. For /i2 = (i.e., ^2 = 0) there is 
no minimum any more. For c? = 4, the disappearance of the minimum above Tc^b for /i2 = 
can be followed in Fig. [61 In the same figure, for h2/hi G {10, 2, 1}, one can see that upon 
decreasing h2/hi the minimum becomes shallower. For c? = 4, the scaling functions d with 
< h2/hi < 2 shown in Fig. [HI exhibit a minimum also below T^^bi which is not seen in the 
data for d = ?>. In (i = 4, below T^^b no minimum can be detected if the surface field hi is 
very strong (see the curves for hi = h2 E {500, 1000} in Fig. [Sj). This suggests, that in = 3 
the surface field hi corresponding to ai = 1 is so strong that surface 1 stays in the regime 
of the normal surface UC (compare the discussion in Subsec. IIII Bp and thus prevents the 
occurrence of a minimum also below Tc^b- Accordingly, a minimum below Tc^b together with 
one above Tc^b is expected to occur for certain values ai < 1 also in d = 3. In order to check 
this expectation, further MC simulation data for these cases are required. 

Also for opposing surface fields the same trends appear in d = 3 and = 4 as evidenced 
by comparing MC simulation data for a film exposed to surface fields corresponding to 



ai = 1 > —02 > 



22j with the present MFT results. 



V. SUMMARY 



We have analyzed the crossover behavior of critical Casimir forces for films of thickness 
L between different surface universality classes. Based on finite-size scaling theory and the 
numerical solution of the corresponding Landau-Ginzburg theory (Eq. ([T])) we have obtained 
the following results: 

1. In Fig. m the scaling function of the critical Casimir force i} {y) is shown as a function 
of the scaling variable y = t (^L/^qY^" with ^ (^^ = ~t~7^ ~^ -'-^) ~ ^o" 
bulk correlation length. In Fig. [H^a) i) (y) is presented for a film with one strongly 
adsorbing surface {hi = 00) and the adsorption preference of the other one continously 
changing from strongly adsorbing the same component of a binary liquid mixture 
(/i2 = 00) to a 'neutral' surface (/i2 = 0) and to strongly adsorbing the other component 
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(/i2 = — oo). The change of the scahng function upon varying the strengths of the 
surface fields for the symmetric {hi = and the antisymmetric [hi = —h2) cases 
is displayed in Fig. dJ^b). A slow convergence towards the strong adsorbing limits 
(+, ±) = {hi = oo, /i2 = ±00) is observed (see also point [3] below). 

2. The properties of films exposed to finite surface fields hi and /i2 can be inferred from 
those of films with (+, ±) fixed point boundary conditions (BCs) by an effective rescal- 
ing (Eqs. f l23p and f l28p ). This scheme is applicable for all films the surfaces of which 
prefer the same component (/ii/i2 > 0). For films with surfaces with opposing prefer- 
ences the scheme is limited to not too asymmetric cases, hi > —h2 > —h2^min > 0, 
where h2,min depends on hi and the scaling variable y. In Fig. [2] the successful corre- 
sponding rescaling of ^ {y; hi, h2) is shown. 

— 1/2 

3. The proposed rescaling predicts the leading order correction ~ hi for the critical 
Casimir force relative to the strong adsorption limit hi = 00 corresponding to (+, ±) 
BCs (see Eq. (!33|) ). The critical Casimir amplitudes vary accordingly as function of 
hi (Fig.E^a)). 



4. Symmetric and weakly adsorbing sur: 



aces, hi = h2 ^ 0, act like an effective bulk field 



/ifo = 2hi/L (see Eq. flM|) and Ref. jl?!). The critical Casimir amplitudes vary ~ —h^^ 
for weak surface fields (Eq. ^S7\ and Fig. El^b)). 

5. The crossover from a purely positive to a purely negative scaling function of the 
critical Casimir force, corresponding to the strong adsorbing limits (+, — ) and (+, +), 
respectively, occurs for the strongly asymmetric case in the sense that the two surfaces 
attract different components of the binary liquid mixture, but one does so much weaker 
than the other, hi ^ —h2 > 0. Within this crossover regime exhibits two maxima, 
one above and one below Tc^b, and the minimum in between is located above Tc^b 
(Fig. H]). The weaker the second surface field is, i.e., for /i2 0^, the weaker and 
broader the maxima are and they are shifted further away from Tc^f,- The amplitude 
of the critical Casimir force in this crossover region is one order of magnitude smaller 
than the maximum value of the scaling function of the critical Casimir force for the 
(+, +) BCs and (+, — ) BCs, respectively. 

In Appendix |B] we have put forward two analytic expressions for the scaling function 
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of the critical Casimir force in the crossover regime in which the rescahng scheme 
(see point |2]) is not apphcable. One (Eq. (IB10|) ) is based on a perturbation theory 
for the order parameter profile around the semi-infinite profile. It reproduces nicely 
the qualitative behavior of in the crossover regime and it is in good quantitative 
agreement for small /i2- The other approach shows that the quadratic interpolation 
of "i? as a function of /12 for fixed y and hi, i.e., d (/12) = a + + ch\, provides 
quantitatively reliable results for values h2.min < h2 < for which 'd (/i2) is not ac- 
cessible by the rescaling scheme. The coefficients a, b, and c are obtained by us- 
ing = [h2) evaluated at three different values (0, h2^s, and h2) of h2 such that 
(/i2 = 0) > -i? (/i2 = /i2,s) = > '/9(/i2 = ^2) where /i2 < h2^min- Figure [31^c) supports 
such an interpolation as it shows that the critical Casimir amplitude A{hi,h2) as 
function of h2 is varying smoothly throughout the whole range of values /i2/^i- The 
performance of both approaches is shown in Fig. |H 

6. For surfaces at which on the one hand ordering is suppressed, e.g., due to missing 
neighboring liquid molecules at the surface, but where on the other hand there is an 
adsorption preference, the interplay of these two opposing influences leads to a richly 
structured {y) (see Figs. |6]and[7]). This structure can be understood by assuming 
that for different values of y the system is characterized effectively by different UCs. 
This suggests that at different temperatures different properties of the surface exert the 
dominant influence on the binary liquid mixture. This is nicely seen for fllms conflned 
by two identical walls which suppress the order (c > 0) and exhibit an adsorption 
preference {hi > 0). If these two influences are comparable, (y) exhibits two minima 
(see Fig. [5]) . The one above Tc^b corresponds to the adsorption preferences and the one 
below Tc^b corresponds to the reduction of the order, as suggested by the comparison 
with the two limiting cases of pure preference and order reduction, respectively. This 
can lead to a situation that the critical Casimir force is attractive above T^^b and turns 
repulsive upon lowering the temperature above Tc^b (see the curve for C2/C1 = 5 in 
Fig. ED- 

7. We have proposed, discussed, and checked a relation between the sign of the critical 
Casimir force and the values of the order parameter at the surfaces relative to its bulk 
value (see Eq. ( HOj) and Fig. |8]). For the asymmetric case (/11/12 < 0) the sign of the 
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critical Casimir force changes above Tc^b if surface 2 is not covered anymore by a layer 
which is enriched in the preferred component. In the case that both surfaces prefer the 
same component {hi > /i2 > 0) and the effect of missing neighbors is not negligible 
(ci,C2 > 0) the critical Casimir force changes its sign twice below Tc^^, approximately 
at the two wetting transition temperatures r^^i and Tw,2 (Eq. (155]) ) of the two surfaces. 
This means that the critical Casimir force becomes repulsive, if at one surface the 
attractive influence of the adsorption preference still dominates and at the second 
one ordering is suppressed relative to the bulk. It switches back to attractive, if at 
both surfaces ordering is suppressed. If one surface is only very weakly adsorbing, 
one change of sign is shifted to a temperature above T^ fe. In the asymmetric case 
{hih2 < 0) the sign changes approximately at the lower wetting temperature, i.e., 
if at both surfaces the ordering is suppressed (see the inset of Fig. [6]). Within the 
present MFT a general, exact, and implicit equation for the temperature at which 
the change of sign occurs is given in Eq. ( !39l) . Accordingly, the change of sign occurs 
at that temperature, at which the semi-infinite profile confined by one surface is not 
disturbed by inserting the other surface. The critical Casimir force is attractive, if 
due to the second surface the ordering is enhanced relative to the semi-infinite profile, 
i.e., if the film enhances the attraction of the component preferred by the semi- infinite 
system. 

8. Similar behaviors of the critical Casimir force as discussed for the spatial dimension 
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- |25| and in Monte Carlo simulation data 



d = A are found in exact results for d = 2 

for d = 3 22| (for d = 3 see Fig. |9]). In all three dimensions, for suitably chosen surface 
properties, the scaling function ■(9 of the critical Casimir force changes sign as a function 
of y. ^ exhibits more than one extrema in the crossover regime and the trends of the 
position of these extrema are the same in all three spatial dimensions. This tells that 
it should be possible to prepare systems in which the sign of the critical Casimir force 
changes forth and back either by changing the temperature or by changing the film 
thickness. It is reasonable to expect that this feature translates to systems consisting 
of colloids the surfaces of which are suitably prepared and which are immersed in 
critical solvents. Therefore the critical Casimir forces provide, at least in principle. 



a mechanism to prepare stable colloidal suspensions. In Ref. [18| it has also been 
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proposed to stabilize colloidal clusters, however by making use of the omnipresent 
quantum Casimir forces. The results we presented here suggest that the use of critical 
Casimir forces provides the potential to control via minute temperature changes the 
distance at which the colloids are in a stable configuration. 
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Appendix A: Short distance approximation for the rescaling functions 

Close to the confining walls the variation of the scaling functions A4(^^,±) {() of the or- 
der parameter profile in a film with fixed point (+, ±) BCs can be well approximated by 



the corresponding scaling function for a single wall in the semi- infinite geometry [37|, i.e., 
|-^(+,±) (C Cw)\ = V^/ (± (C ~ Cw)), where E {±0.5} denotes the position of the sur- 
face, ± (C — (w) corresponds to ( > = —0.5 or ^ < = 0.5, and the sign of Ai equals 
the one of the corresponding (infinite) surface field. If a surface field hi, i G {1, 2}, is strong 
enough this relation can be used for that surface i in Eqs. ( l25i) and ( 126|) which determine 
the rescaling functions r and Co- If both surface fields are strong, one obtains 

rid — 2 = 1^1 1 (Al) 

(r.4-0.5-a'^)+0.5)^ 

/2 

r^d 2 = \h2\- (A2) 

^'(-r.40.5-a'^)+0.5)' 

Equations (jAl]) and (jM]) can be written as 2^/^/ \hi\^^^ = -0.5 =F Co"' + 0.5/r,d, with =F - 
for i = 1 and =F — !■ + for i = 2, from which one obtains Eqs. (13T!) and (!32l) . Within 
this approximation the contributions of the two surfaces to the rescaling parameters are 
independent of each other. Moreover, within this approximation Eqs. (pTl) and (1321) offer 
a transparent interpretation of these functions: Co is a measure of the asymmetry of the 
surface fields the film is exposed to, and r provides a measure of the deviation from the 
fixpoint BCs (+,±). 
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Appendix B: Analytic expressions for (y) in the crossover regime 

We propose two analytic expressions for the scaling function of the critical Casimir force 
in the crossover regime, i.e., in which the rescaling scheme is not applicable. One is based 
on a perturbation theory of the order parameter profile around the semi-infinite profile. The 
other approach is the quadratic interpolation of "t? as a function of h2 for fixed y and hi. 



1. Perturbation theory for the critical Casimir force 

In the crossover regime the scaling function of the critical Casimir force exhibits a 
rich structure with the emergence of three extrema and up to two changes of sign. In 
order to capture these features qualitatively and quantitatively by an analytic expression, 
we approximate the actual order parameter profile $ by a term $o which satisfies the ELE 
(Eq. (jl])), but in general not the BCs in Eqs. ([5]) and IQ, and a perturbation part 5$: 

$ (z) = $0 (z) + 5$ (z) . (Bl) 

Inserting the ansatz of Eq. flBl|) into Eq. (|4]) provides the differential equation determining 
the deviation 5$ = $ — $0, i-e., 

= (^r + |<l>2) 5$ + (9 ((5$2) , (B2) 

where 5$' = 85^ /dz and O (5$^) stands for terms quadratic and cubic in 5$. The BCs 
are (compare Eqs. ([5]) and ([6]), for simplicity we consider here and throughout this appendix 
only the crossover from the special to the normal transition, i.e., Cj = 0) 

= -V-'^o(-V2) = -AV (B3) 

L/2 



dz 
and 



= ho,2 - $0 iL/2) = Aho,2. (B4) 

z=L/2 



dz 

(If $0 happens to be the full solution, 5$ = and thus A/io,j = 0.) The critical Casimir force 
is then obtained from the ansatz $0 + (with 5$ (r, /io,i, [$o]) determined by Eqs. (lB2p - 
(iBiD ) via the stress tensor (Eqs. ([7]) and ([8])): 

fc = fc (r, [%]) = %z [$o + m - %z I'^b] 

= %z [$o] + $[,5$' - (r$o + |$o) 5$ - %z [$5] + O , (B5) 
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i.e., fc is a function of r and a functional of $o- Here and in the following we do not 
indicate explicitly the dependence on the surface fields /io,i- Concerning the critical Casimir 
force Oj^the critical temperature, in the following we shall neglect the explicit dependence 
of fc on r and assume that it enters only via the profile $0) i-e-? fci'^ 7^ 0, ['I'ol''")]) ~ 
/c(r = 0, [% (r)]). 

At the temperature Tc, at which the critical Casimir force is zero and changes sign, the 
actual profile $ in the film {—L/2 < z < L/2) coincides, in the case /lo.i > I ^0,2! assumed here 
and in the following, with the corresponding semi-infinite profile $* (2; > — L/2, Tc) because 
at this temperature there is no finite-size contribution to the free energy (see Subsec. IIII Cp . 
We are interested in the change of sign and in the crossover regime around Tc where the 
critical Casimir force is small. For this case and by invoking the above approximation 
scheme it is reasonable to dispose of the not yet specified $0 such that $n(r) = $s (r). 



where for the time being $5 (r) = $s {z, r; 2:0) is any (analytically known [4^) semi-infinite 
solution of the ELE; zq = zq (/io,i, z^) is a lengthscale which encodes the dependences on the 
surface field /io,i and the position z^ of the confining wall. The specification of this zq will 
be discussed below. 

In order to proceed we now determine 5$ explicitly by assuming that it is sufficiently small 
and, accordingly, we neglect in the following terms O (5$^). For reasons of simplicity we re- 
strict ourselves to solving the linearized ELE (Eq. flB2|) ) at bulk criticality. According to the 
above approximation scheme one has $0 (^; r = 0) = $s (2;, r = 0; Zq) = y^l2/g {z + -2o)~\ 
so that the linearized Eq. (lB2p turns into (5<l>" {z) = 6{z + zo)^'^6^ with the solution 

5$ [z; T = 0)=Bi{z + zof + B2{z + ^o)"' • (B6) 

The coefficients Bi and B2 are determined by the BCs in Eqs. flB3l) and flB4l) : 

5i = i {z~^Aho,i + zZ'Aho,2) B (B7) 

and 

^2 = ^ (4A/io,i + zlAho,2) B, (B8) 

with B = z\z^_ {z\ — 2;^) ^ and z± = zq± L/2. 

The value of the stress tensor for $0 = equals its bulk value, i.e., Tzz [^s\ —Tzz = 0, 
and within the linearized theory (i.e., neglecting terms O (5$^)) the critical Casimir force 
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at r = is approximated by (see Eq. ( IBSP ) 



/c(r = 0,[<l>,]) = -5v/l27^fii. 



(B9) 



Straightforward calculations show that fc does not depend on 

We now resume the remaining task to specify the value of the parameter zq. We have 
chosen the straightforward and simplest condition A/io,i = 0, i.e., we consider for $o = the 
semi-infinite profile corresponding to surface 1. The condition A/io,2 = is not considered, 
because we have observed that the former choice yields better results and it ensures that <l>s 
is non-singular throughout the film for all finite values of hi, i = 1,2. Within the crossover 
regime the profiles $3 + obtained in such a way compare well with the ones obtained by 
the full, numerical minimization of H, while the critical Casimir amplitude A as function of 
/i2 can be approximated in this way only for values /i2 ~ /i2,s (where A = for /i2 = /i2,s)- 
(The full MFT profile $ minimizes the Hamiltonian Ti. Therefore a more sophisticated 
way to determine zq would be to treat it as a variational parameter and taking that value 
Zq = ^o,m which minimizes ($s {z; Zq) + 5$ (z; Zq)). ) 

In order to extend Eq. (1B9I) to values r 7^ 0, we apply the approximation stated below 
Eq. (]B5|) . i.e., we take into account only that dependence on r which enters via the profile $o- 
To this end we express the right hand side of Eq. (IB9|) in terms of the profile $s [z] r = 0). 
This can be done by using the relations $5,(1,2) (r = 0) = 12 /gz~^ and ^-^ 2) ('^ = 0) = 
— ^/12/gz^'^, where = ^s{z = zq^ L/2; r = 0) is the value of the profile $0 = at 
the surface i and =F ~ for i = 1 and =F ~^ + for i = 2. Expressing Bi in terms of ^s,i 
and j is not unique, yet it is restricted by the condition to preserve the relation fc ~ 
(see Eq. f|T5l) and the note after Eq. fl20|) ) and to preserve the symmetry with respect to 
interchanging 1 and 2. 

We have found empirically that the final result is not sensitive to the particular way of 
expressing Bi in terms of ^s,i and ^ and in the following we choose a replacement in which 

j replaces only the terms z"^ in the denominator of i? = 1/ [^z'^zZ^ — z'^z^'^Y Together 
with Eqs. flB3l) and flB4p and in terms of the scaling variables given by Eqs. f|TO|) - f[T^ and 
(!20|) this choice leads to the approximation 



Ms,i = Ms {(= -0.5; y) and Ms,2 = Ms{C = 0.5;y) [M',^^ = f^Ms {( = -0.5;y) and 
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^{y;hi,h2) 



99 Ml,{M's^,)-'-Ml,{M's^,)-' 



(BIO) 



^ Ais (C = 0-5; y)] are the values at the two confining surfaces of [the derivative of] 



'■,2 ~ dC 

the scahng function corresponding to a semi-infinite order parameter profile: 

-1 



[sinh (yV2 (^^^^ 



c + c 



i/>o 

y<0 



(Bll) 



with the parameter C, corresponding to Zq (see above). We have used for C, that value C, {y, hi) 
which follows from the condition Ai'^ i = —hi (i.e., the one corresponding to that value of zq 
as chosen above). From Fig. H] we infer that this approximation (Eq. flBlOp ) captures nicely 
the qualitatively behavior of i} (y) in the crossover regime where three extrema emerge and 
changes sign. 

Within this regime, for small \h2\ it is also in good quantitative agreement with 
d {y\hi,h2) as obtained from the full numerical minimization of "H (see the curves for 
^2/^1 € {0, —0.003, —0.007, —0.01}), whereas for stronger /i2 it underestimates 'd. From 
data not shown we infer that the range of values of /i2/^i, for which the expression in 
Eq. flBlOP (with Q chosen as above) approximates well (y; hi, h2), increases with decreas- 
ing hi, e.g., for hi = 0.25 it captures {} even for /12 = ±hi. A more sophisticated choice 
for which takes into account not only the value of hi but also the one of h2, may even 
improve the value for given by Eq. (IBlOp and may extend the range of applicability of 
this approximation. 



2. Interpolation of the scaling function in the crossover regime 

As discussed in the main text, apart from the crossover regime > /12 > /i2,mm (see 
Eq. (!29|) and the discussion before it, considering here and in the following hi > \h2\ > 0), 
the scaling function [y; hi, /i2) has the same functional form as ■})(+ +) (see Eq. fl28l) ) and 
its value is known analytically, yet implicitly, because both '&{+,±) [3(| and the rescaling 
function r = r {y; hi, h2) (Eqs. fl2^ and fl2Bl) ) are known. In Appendix IB II we have provided 
a perturbation theory for the scaling function of the critical Casimir force. It results in an 
explicit expression for ^ but it is, depending on the value of hi, limited to a certain range 
of values of /i2 and does not necessarily apply for the whole crossover regime. 
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In the following we present a calculation scheme for ■(9 which covers the whole crossover 
regime and which avoids the full numerical minimization of % (Eq. ([1])). It turns out that 
a quadratic interpolation of 'd as function of h2 for fixed hi and i.e., 

d'''°''{y]hi,h2) = ahl + hh2 + c (B12) 

adequately serves this purpose. The coefficients a, 6, and c which depend on y and hi are 
determined by the condition that ■O'^'^"'^'^ takes the (exact) values of d at three distinct values 
of /i2 5 for which d is known in terms of analytic expressions, and which we shall specify in the 
next step. In order to interpolate the value of d in the crossover regime we take one value in 
each region where the rescaling scheme applies. As the first value we choose /i2,+ = 0, because 
it is at the boundary of the crossover regime and it simplifies the calculation of d (as given 
by Eq. (1251) ) because for h2 = Q Eq. fl2Bl) immediately yields Co = 0.5 and one is left with only 
one implicit equation, i.e., Eq. ( l25l) . As the second value we choose h2- = h2^min {v = 0) 
for all values of y. This choice avoids that one has to calculate the corresponding h2,min (y) 
for each value of y. (For h2- = h2^min {y 7^ 0) the ensuing interpolation intervall would 
be narrower than for other choices of /i2,- and thus one would expect in that case more 
accurate results for {}^'''°'^'^, However it turns out that h2^min {y = 0) is sufficiently close to 
h2,min {y 7^ 0) to yield satisfactory results; we recall that h2^min (?/ = 0) < h2^min (?/)•) As the 
third value we take /i2,s = M'g (C = 0.5; y, hi) < for which the critical Casimir force is zero 
and thus ^ (/i2,s) = (see Eq. fl39l) ). Since the rescaling procedure amounts to stretching 
the fixed point scaling function ±), it cannot describe the qualitatively different shapes 
of i) in the crossover regime, such as the change of sign of i}. Accordingly, h2^s lies within 
the crossover regime. In sum, the conditions fixing a, b, and c are: 

r™^^ {y; hi, h2 = h2,+) = [r (y; h, /i2,+)]' ^9+,+ ([r (y; hi, /i2,+)]"' v) , (B13a) 
^9--^ (y; hi,h2 = h2,.) = [r (y; h, h2,^)f ([r (y; h, /12,-)]"' v) , (B13b) 
r™^^(l/;/ii,/i2 = M = 0, (B13c) 

with /i2,± and /i2,s as given above and r = r {y; hi, h2) follows from Eqs. ( l25i) and ( l26l) . The 
performance of (Eqs. (1B12I) and (1B13I) ) is in good agreement with the data obtained 

by the full, numerical minimization of "H (Eq. (Q) (see Fig. Hj). It is worth mentioning, that 
a linear interpolation (i.e., a = 0), using only two of the conditions in Eq. flB13p . performs 
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poorly. 
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FIG. 1: The scaling functions ?? of the critical Casimir force are shown for various types of confining 
surfaces covering the crossover between the normal transition (it, i.e., hi = ±00, strong adsorption) 
and the special transition {sp, i.e., hi = 0, neutral surface, ci = C2 = 0). In (a) /ii = 00 is fixed and 
/i2 varies from +00 [(+,+)] to —00 [(+,—)]. For all cases shown, t9 for finite /i2 can be expressed 
in terms of i9(+,±) by using Eq. (j28p . In (b) is shown for the symmetric and antisymmetric cases 
hi = I/12I and —00 < /12 < 00. In general the convergence towards the strong adsorption limit 
(black lines) is rather slow (compare Fig. [3]^a)). 
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FIG. 2: For all surface fields with /11/12 > and certain regimes of surface fields with /11/12 < (see 
Eq. (j29p and the discussion before it) the scaling function i) of the critical Casimir force for finite 
surface fields can be mapped onto the fixed point scaling functions and _), respectively: 

r"^-!? (y; /ii, /12) = ±) ^ = ^ {y, hi, (Eq. ([28]) ). The symbols correspond to hi = 168 

for different values of /i2 obtained by minimizing the Hamiltonian in Eq. ([TJ and r (y; /ii, /i2) is 
given implicitly by Eqs. (j25p and (j26p . and are known analytically jsol] and given by 

the full black lines. Here d = 4 and = 1/2. 
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FIG. 3: In (a)-(c), symbols represent data obtained from the numerical minimization of the Hamil- 
tonian in Eq. ([T]) and the full lines follow from the rescaling scheme (Eq. (j28l) ). (a) For /ii — t- oo 
the critical Casimir amplitude A-(- (hi) = A {hi = /12) = (y = 0; /ii, /i2) approaches its fixed point 
value _|_) slowly ~ h^ from above (see dotted line representing — (l/2)x -|- 1.1 and Eq. ([33]) ). 
(b) For small surface fields hi = /i2 and ci = C2 = the critical Casimir amplitude A_|- (hi) 
varies ~ —h^ (see the dotted line representing (4/3) x — 0.8 and Eq. ()37p ). (c) The variation of 
A (/ii, /12) as a function of /12 for fixed hi € {1, 168}. The values are normalized by |A+ {hi)\ with 
A+ {hi = 168) / |A(+^+)| ~ -0.51 and A+ (hi = 1) / |A(+^+)| ~ -0.010. Figure EJc) focuses on 
small values of h^ for which the Casimir amplitude changes sign and for which the rescaling prop- 
erty (Eq. (j28p ) does not hold. In this latter regime the full lines do not apply and are replaced by the 
numerical obtained data points. These data demonstrate that also in this regime A (/ii, /i2) varies 
smoothly. For /12//11 1 both curves approach —1 by construction and due to the fact that for 
hi = /1.2 the Casimir amplitudes are negative. Upon further increasing /12 the Casimir amplitudes 
turn more negative and approach the limiting values A (hi = 168, /12 = 00) / |A(_|_ _|_-)| ~ —0.70 and 
A (hi = l,/i2 = 00)/ |A(+_+)| ~ -0.099. For hi = 168 and /loAi -1 the full line gets close 
to the MET value A(_|_ _)/ |A(_|_ _,_)| = 4 known for \hi\ = 00 30(]. This property holds also for 
hi = 168 in the limit /i2 — —00, i.e., A{hi = 168, /12 = —00) /\A{hi = 168, /i2 = cxd)| ~ 4. For 
hi = 1, A accidently has the same strength for /12//11 = ±1- For more negative values of /i2 < —hi, 
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A {hi = l,/i2) first increases. However, upon further decreasing /12, it will decrease 
and change sign at /i2 = —39.5 (as obtained from Eq. (j39p ) and remains negative with 
a slightly negative limiting value for /12 — >• —00. This can be infered from the fact that 
A(/ii,/i2 = —00) as function of hi changes sign at hi^s = (Eqs. (|39]) and (|Blip ) 
and is negative for hi < hi^s- 
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FIG. 4: The scaling function ?? of the critical Casimir force changes from being purely positive to 
being purely negative by the formation of three extrema in the crossover regime for ci = C2 = 0, 
fixed /ii > {hi = 168 in the figure), and /12 — t- 0^. Symbols correspond to data obtained by a full, 
numerical minimization of the Hamiltonian in Eq. ([1]). For hi = 168 the values /i2/^i = 0, —0.025 
represent the limits of the range of values of /12 for which the rescaling scheme (Eq. (j28p ) holds 
(i.e., 112/ hi ^ (—0.025, 0)). The qualitative behavior of ?? for h2/hi G [—0.025, 0] is captured by the 
description obtained from the ELE linearized in terms of the deviation (5<I> = <I> — <l>s of the order 
parameter profile <I> from the profile in a semi- infinite system (see Eq. (jBlOp . full lines). For small 
/i2 it is even in good quantitative agreement (see curves for h2/hi G {0,-0.003,-0.007,-0.01}). 
Values obtained by the quadratic interpolation scheme given in Eq. ([BT2]) provide a very good 
description of the scaling function in the whole crossover regime (dotted lines for h2/hi G 
{-0.025,-0.02,-0.015}). 
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FIG. 5: The scaling function i!) of the critical Casimir force for fixed ci = C2 = 100 and in- 
creasing surface fields hi = h2 crosses over from the behavior typical for the {ord, ord) BCs 
(ci = C2 = oo, /ii = /i2 = 0) to the behavior of the (+, +) BCs. Thereby the minimum below T^^ 
disappears and a minimum above Tc^f, is formed. For moderate values of hi there are two min- 
ima (see, e.g., hi = 168). For comparison the scaling function 'd[ord.ord) {u) (black line) is also 
shown. ??(+,+) {y) is off the scale of the figure ('i?(_(_^_|_) {y = 0) / |A(_|_ _|_)| = 1). Notice that, even for 
hi = 1000, 'd {y; ci = C2 = 100, /ii = /12 = 1000) differs markedly from This is probably due 

to the strong surface enhancement parameters q and the slow convergence apparent in Fig. [1] 
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FIG. 6: For fixed surface parameters (ci = 100, hi = 168) of one surface varying the surface field 
/i2 of the second surface with fixed C2 = ci induces a rich variation of the scaling function -d of the 
critical Casimir force. Whereas for /i2/^i ^ —4 the characteristics of resemble those of ■!?(_!_ _), 
for /i2//ii = 10 the scaling function exhibits the characteristics of and ^(^^ord) above and 

below Tc^b, respectively. For — 3 < /i2//i2 < 5, exhibits an even more complicated structure (see 
the main text, Subsec. flllBp . The range of values for y around those corresponding to the wetting 
transition temperatures, y^,^! = —2.8 for surface 1 and yw,2 = —11.3 for surface 2 with \h2/hi\ = 2, 
is shown enlarged in the inset for the curves belonging to /12//11 G {-2,-0.9,0,2}. As discussed 
in detail in the main text (Subsecs. IIIIBI and IIIICp . the change of sign occurs up to a finite-size 
correction 5y at the wetting transition temperature y^ if the film crosses at this temperature from 
a UC with an attractive force to a UC with a repulsive force (or the other way round). For 
^2/^1 S {~2,0, 2} the correction is rather small, 6y < 1, but 6y = 2.6 for /i2/^i = —0.9. 
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FIG. 7: The scaling function -i? of the critical Casimir force for the case that one surface is fixed 
at (ci = 100, hi = 168) and at the other surface the same surface field is applied, /12 = /ii, but the 
surface enhancement C2 is varied. The structure in the variation of "i? (y) can be understood by 
assuming that the film crosses asymptotic regimes of different fixed points upon varying the scaling 
variable y given in Eq. (I12p (see the main text). As discussed in the main text (SubsecUlB]), the 
behavior is similar to the one shown in Fig. [6j The case C2/C1 = 5 shows that it is possible that 
above T^.^ the critical Casimir force is attractive and that upon lowering the temperature it turns 
repulsive above Tc^b. 





C2/C1 = 50 : 


. . . 1 . . . 


l\/o.5 /y 

. . . 1 ... 1 ' 



42 




1 — ' — ' — ' — I — ' — ■ — ' — I — ' — ' — ■ — I — ' — ' — ' — r 




-40 -20 20 40 



FIG. 8: Analysis of the changes of sign of the scahng function of the critical Casimir force as a 
function of the scaling variable y = tL'^ (Eq. (jl2p ) for two examples. Figure [H^a) shows a typical 
example for the crossover regime between the normal and the special transition. The values of the 
surface parameters are ci = C2 = and hi = — IO/12 = 2.5. Figure [5]^b) represents the crossover 
regime between the normal and the ordinary transition. Here the values of the surface parameters 
are ci = C2 = 100 and hi = IO/12 = 1680. In the insets the scaling functions Ai (C, y) of the order 
parameter profile in the film are shown for those values of y indicated in the main figures. The 
corresponding bulk values Aib{y) (Eq. p^ ). which are for y > 0, are indicated as horizontal 
dotted lines. One can infer from these plots that changes sign approximately at that value of 
y at which the value of the order parameter M. {(, y) at one surface equals its bulk value M.^ (y) 
(compare Eq. (00])). 43 
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^'IG. 9: The scaling function -i? of the critical Casimir force obtained from Monte Carlo simulations 
22l | for a d = 3 Ising film of size >i Ly x L^, with Lx/3 = Ly/3 = = 10 and with periodic 
BCs in the x and y directions. The surface field hi of surface 1 is fixed and corresponds to ai = 1 
(see Subsec. flV Bp and the surface field /12 of the surface 2, corresponding to 02, is varied. Upon 
decreasing /i2 the minimum becomes shallower and a maximum emerges below T^.^ (see 02 = 0; 
for 02 = 0.1 there exists a very weak and broad maximum which is hardly seen on the scale of 
the figure). This behavior is also observed for d = 4 with surfaces in the crossover regime between 
the ordinary and the normal transition (see Fig. [6] and the main text). Here y = t {L/f^)^^" with 



L = L,,t = {T- T,,,) /Te,,, ^ 0.50 



461], 1/ ~ 0.63 [45], and A(_^^_^) (d = 3) ~ 0.38 
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